** Mark D. Ramirez and Reed M. Wood
** Do File 
** Authoritarian Opposition? Authoritarian Disposition and Resistance to Public Health Mitigation Strategies during COVID-19
** Political Research Quarterly (Forthcoming) 


** Set WD
use "PRQ_AUTH_USUK.dta"

set more off

svyset _n [pweight=weight], vce(linearized) singleunit(missing)


** Baseline model US
svy: reg zPROTEST authoritarian_score trust use_soc_media vax_belief_harm Partisanship Ideology Education age Gender rel_id urban if US==1

** State fixed effects (Figure 1; Panel 1; Model 1)
svy: reg zPROTEST authoritarian_score trust use_soc_media vax_belief_harm Partisanship Ideology Education age Gender rel_id urban ib12.inputstate if US==1

estimates store m1 

margin, dydx(*)

** Baseline model UK
svy: reg zPROTEST authoritarian_score trust use_soc_media vax_belief_harm Partisanship Ideology Education age Gender rel_id urban if US==0

** Region fixed effects (Figure 1; Panel 1; Model 2)
svy: reg zPROTEST authoritarian_score trust use_soc_media vax_belief_harm Partisanship Ideology Education age Gender rel_id urban ib7.profile_gor if US==0

estimates store m3

margin, dydx(*)

** Figure 1; Panel 1
set scheme s1mono
coefplot (m1, label(US) pstyle(p1)) ///
		 (m3, label(UK) pstyle(p2)) ///
		 , keep(authoritarian_score trust use_soc_media vax_belief_harm Partisanship Ideology Education age Gender rel_id urban) xline(0) title("Protest Support")  

		 
** Coefficient test		 
suest m1 m3

test [m1]authoritarian_score=[m3]authoritarian_score		 
		 
	
** US ideology interaction 
svy: reg zPROTEST c.authoritarian_score##i.Ideology trust use_soc_media vax_belief_harm Partisanship Education age Gender rel_id urban i.inputstate if US==1
estimates store m2
margin, dydx(*)
margins Ideology, at(authoritarian=(0 .25 .5 .75 1)) vce(unconditional) level(90)
marginsplot
		 
** UK ideology interaction 
svy: reg zPROTEST c.authoritarian_score##i.Ideology trust use_soc_media vax_belief_harm Partisanship Education age Gender rel_id urban i.profile_gor if US==0
estimates store m4
margin, dydx(*)
margins Ideology, at(authoritarian=(0 .25 .5 .75 1)) vce(unconditional) level(90)
marginsplot


gr combine US_protest.gph UK_protest.gph, iscale(.5) ycommon



